General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA 

Technical Memorandum 78090 


Theory of Flux Anisotropies 
in a Guiding Center Plasma 

(NASA-TH-78090) THECET CF FIUX ANISCTECPIES N78-21028 

IN A GOIEING CENTEF P1ASNA (NASA) 19 p HC 
A02/NF A01 CSCl 03E 

Unclas 

G3/90 12572 


T. J. Birmingham and T. G. Northrop 


FEBRUARY 1978 


National Aeronautics and 
Space Administration 

Goddard Space Right Center 

Greenbelt, Maryland 20771 




TM 78090 


THEORY OF FLUX ANISOTROPIES IN A GUIDING CENTER PLASMA 


T. J. Birmingham 

Laboratory for Extraterrestrial Physics 
and 

T. G. Northrop 

Laboratory for High Energy Astrophysics 
Goddard Space Flight Center 
Greenbelt, Maryland 20771 


February 1978 


GODDARD SPACE FLIGHT CENTER 
Greenbelt, Maryland 


THEORY OF FLUX ANISOTROPIES IN A GUIDING CENTER PLASMA 


T. J. Birmingham 

Laboratory for Extraterrestrial Physics 
and 

T, G. Northrop 

Laboratory for High Energy Astrophysics 
ABSTRACT 

Assuming time stationarity of the one particle distribution function f 
on the scale of the bounce motion of particles in a magnetic field B, we 
expand the Vlasov Equation through 0(t) in the adiabatic parameter e, 
which is the ratio of particle gyroradius to scale length of the magnetic 
field. Since f is directly proportional to particle flux dd>/dWdf2 differential 
in kinetic energy W and solid anple f is in principle measurable in space 
experiments, and our analysis is tailored to be explicitly applicable to 
space problems. To 0(1), f is gyrotropic; its first velocity moment is (if 
non-vanishing) parallel to B, and hence macroscopic parallel flow is 
included in this term. The 0(e) contribution is non-gyrotropic and 
macroscopic flow 1 to B plus additional parallel flow results from these 
terms. The degree of non-gyrotropy and hence the amount of cross-field 
macroscopic How depend on the perpendicular component of the electric 
field E, on curvature and shear in the magnetic field, and on the spatial 
gradient f Q , pitch angle derivative df 0 /35 and speed derivative 9f 0 /dv of 
the lowest order distribution function f Q . We also show that the usual 
expression for the electric field E wliieh produces plasma co-rotation in an 
axisymmetric system such as a dipole also holds for any non-axisymmetric 
but rigidly rotating magnetic field pattern, provided the observed magnetic 
field is used in place of the dipole field. 


ui 


f BDLK&IBG PAGE BLANK. NOT fl UM 

CONTENTS 

Page 

ABSTRACT ili 

I. INTRODUCTION 1 

0. EXPANSION OF VLASOV EQUATION 3 

ni. THE ELECTRIC FIELD IN A RIGIDLY 

COROTATING MAGNETOSPHERE 13 

REFERENCES 15 


v 


THEORY OF FLUX ANISOTROPIES IN A GUIDING CENTER PLASMA 
I. INTRODUCTION 

We often consider charged particles interacting with the electric and magnetic fields 
E and B in a magnetospheric environment as behaving adiabatically: time variations are 
sufficiently slow and spatial variations sufficiently long that individual particles undergo 
nearly periodic motions. Under adiabatic conditions the Vlasov Equation 
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which governs the distribution function f(£, v, t) of such particles, can be expanded 
asymptotically in the adiabatic parameter. 

In this paper (Section II) we carry out a formal adiabatic expansion of the Vlasov 
Equation in a manner paralleling that of Hastie et al. [1967J, but we tailor it to be mors 
explicitly applicable to space problems, such as the one in the following paper. We do so 
assuming that f is stationary on the bounce time scale of particles moving in B. We fmd 
that f is constrained in form: the dependence of f on X, the gyrophase angle of v with 
respect to B, is related to the spatial gradient vf Q of the lowest order f, to the deriva- 
tives 9f 0 /3v and df 0 /85 with respect to speed and pitch angle respectively, and to the 
magnitude, curvature, and shear of B as well as the magnitude and direction of E. This 
X-dependence of f is called gyro-anisotropy and leads, for example, to a net particle flow 
perpendicular to B. 

Given values of E and B and of the derivatives of f 0 , our result can be used to cal- 
culate f and its velocity space moments for comparison with charged particle measure- 
ments. Measurement of the flux of particles dd>/dWdS2 differential in energy and solid 
angle is equivalent to measuring f because of the well known relationship 
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dWd« m 2 (2 > 


1 


Conversely we may deduce the values of one or more of the sources of gyro-anisotropy, 
i.e., of the X-dependence, by examining the phase space variation of measured fluxes. In 
a subsequent paper, to which this paper is a formal theoretical prelude, we apply this latter 
approach to Pioneer measurements of energetic protons in Jupiter’s magnetosphere. 

The electric field E v hich enforces plasma “eorotation” is the most generally agreed 
upon E to which majpictospheric particles are subjected. In Section 111 we derive for use 
in our f equation the functional dependence on r and t of the corotational E for the 
case of a magnetic rotor whose magnetic field has no symmetry axis. We require that the 
magnetic field pattern, whatever it might be, rotate rigidly, so that the field at any time 
looks like that at an earlier time, but rotated. Our work here is an extension of previous 
theoretical work of Mestel [1961, 19681 . 


II. EXPANSION OF VLASOV EQUATION 


' Let us define a triad of locally orthogonal unit vectors 6 j(£, t), e 2 (£, t), and e 3 (£, t) 
with !| s B/|B| and e 2 and oriented arbitrarily (in this paper) in the plane transverse 
to B. Because we usually measure particle fluxes as functions of the energy and of the 
direction of the velocity vector, it is convenient to introduce the particle pitch angle 


and its gyrophase angle 





(3) 


( 4 ) 


and to write the Vlasov Equation using v, 6 , and X as variables. 
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In Equation (5), E } , E 2 , and E 3 are the components of E in the directions of our unit 
vectors. Note that the spatial gradient yf in (5) is at constant (v, 6 , X, t) and differs from 
that at constant v because the reference directions (e t , § 2 , and ), with respect to which 
5 and X are defined, change with £. However, f, even though expressed as a function of 
v, 8 , X, and £ from here on, is still the density in £, v space. 
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In the adiabatic expression we write 


f - f 0 + f t + • • * 


where successive terms become smaller by the ratio 
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of the particle gyroradius r Q = v/(q| B|/mc) - v/w to the characteristic scale length L of 
the magnetic field. Note that € depends on particle energy as W h . In Jupiter’s dayside 
mapwto^nhere between 20 and dO Jovian radii from the planet e is typically ^.01 for a 
1 MeV proton. There are magnetospheric regions, however, where B is very weak as in 
the Earth's tail near the neutral sheet or L is very small as in the current sheet region in 
Jupiter’s dawn meridian (Smith et al. 1 19761, Goertz et ai. (19761) and adiabatic theory 
is invalid for a 1 MeV proton. 

We make the small electric field assumption that E - 0(e) for both the parallel Ej 
and perpendicular components E 2 and £ 3 , The smallness of E, is necessary in order that 
particles do not freely accelerate along B, thereby destroying any semblance of periodicity 
and adiabaticity. Whether E 2 and E 3 are small of 0(e) depends on the ratio of the elec- 
tric field drift speed cE^/B to the particle speed v. In the following paper we take E 1 to 
be the electric field which produces corotation and find that this ratio is about 10" 3 p for a 
1 MeV proton, where p is the distance from Jupiter’s spin axis in Jovian radii. Since p 
<50, cE A /Bv does not exceed .05, and hence it is indeed permissible to choose E of 0(e) 
for the purposes of the following paper. 

The relative asymptotic sizes of terms in (5) can now be estimated. A handy device 
here is to regard the coefficient m/q as carrying an intrinsic weight of e. All terms in (5) 
are at largest of 0(1) except the first and last, and to lowest (1/e) order 


Equation (6) expresses the result that any gyrophase dependence of f 0 gyrates about B 
at the gyro-frequency to. We demand stationarity on the gyro-timescale so that 
df 0 /3t * CKO, 3f 0 /dX s 0, and f 3 = f 0 (v, 6, r, t). 

We now demand further that 3f/3t s 0(e)» Le., that f be stationary on the time scale 
of the bounce motion. With this restriction there ean be no particle bunches moving 
along B, The 0(1) contribution to (5) thus becomes 
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which may be integrated to yield 
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where v 1 = v / v. The integration in (8) can be carried out by noting 


v s e. cos 5 + sin 6 (e, cos X* + e^ sin X'i 
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The result is 
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Since X it an angle variable, f must be periodic in it f 0 , being independent or X, is 
dearly so. Everything in Equation (9) for f 3 is periodic in X with the exception of the 
Unear terms. We enforce periodicity by constraining f 0 further so that the coefficient of 
X vanish®: 
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Now note that 

0 ® V • B s V * (f , B) s e, • 7B + [(e, * Wj ) * # 2 + (e 3 • V4j ) * e 3 ) B (11) 
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Let us denote e, • V by 3/3s. Equation (12), being linear in f Q , is satisfied by any 
arbitrary function of the characteristic constants which label flow lines in the 3-dimensional 
s, 5. v space. These constants are obtained by solving the ordinary differential equations 
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From the flrst pair of these equations comes the parallel energy 


as one constant, and from the second pur we find (uang 14) the magnetic moment 
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to be a second. The constraint (10) therefore dictates that on any given field line 
f 0 = f 0 (K, M). Because M is double-valued in 5, separate distributions may be specified 
for particles moving parallel and antiparallel to B (Hastie et al. 1 1967] ). If the particles 
are trapped in a magnetic mirror, the stationarity of f requires that these two distributions 
be the same. However, there are situations, sueh«as circulating particles in a toroidal de- 
vice or particles precipitating at a steady rate in a magnetosphere, where the symmetry 
need not exist. In addition, the form of the dependence of t' 0 on K and M may differ on 
different field lines. If we label a field line by the parameters a and p (Northrop 119611. 
Stem 119701) then ( Q = f 0 (K. M, a. $). 

Note further that the remaining Vin.depc ndent terms in (9) can be lumped together 
by redefh ing the arbitrary function f, . Thus, through CKe), 

f ® f 0 (K. M.a.p' + f,(v. 5. r) 
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Equation (16) is the principu result of this paper. Given f Q + f, and the magnetic 
and electric fields B and E, it a. lows one to calculate the non*gyrotropy (which contrib- 
utes, for example, to the cross-field particle flow). The existence of f t means that the 
gyrvV«*oic part of f may have a small CHe) component which depends in an arbitrary 
fed- to* on s, v. and 5 and still maintain df/dt = 0 + 0(e). If the time variation of f were 
limited yet further, restrictions on f, . similar to those imposed on f Q here, would arise 
(Hastie et al. (1967)). It should also be noted that because of the redefinition of f r 
f 0 + f, is no loneei f(X = 0) but is the X average of f. 

It is instructive to examine velocity space moments of (16) to see the descendants of 
each term. The density 

n ~ K J d 3 '■*' d J v(f 0 + f, ) * / d 3 vf 0 + Q(<r) (17) 


does not involve any X-dependent terms, because they integrate to zero. To the macro- 
scopic flow velocity 
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on the other hand, f 0 + f, contributes the parallel (to B) component and the cos X and 
sin X terms in f are responsible for the perpendicular fluxes in the 4 2 and directions 
respectively. In order that a non-zero < v > * 4, exist, (f 0 + f, ) for particles with pitch 
angles 0 < 6 < *12 must be different from (f G + f, ) for particles with pitch angles 
x/2 < 6 < tr. 


An experimental detector frequently samples only a portion of the An steradians 
through which v may vary. For example, in the Pioneer experiment analyzed in the 
following paper, the proton detector scans in 6 through 2 n radians in a plane to which 
B is arbitrarily inclined. An experimental anisotropy is then obtained from the magnitude 
and direction of the quantity 



dfl ftv, 6<0),X(0)1 y(0>. 


The cos 2X and sin 2X terms in (16), which vanish upon integration over ail angles of v, 
contribute lo this “reduced How velocity” provided the pitch angle distribution of par- 
ticles is non-uniform so that 3f 0 /3S # 0. 

The perpendicular flow terms in (18) can be cast into a more familiar form by some 
straightforward integrations, integrations by parts, and vector manipulation. Details are 
omitted because we are interested only in the fate of terms in (18). 
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+ X (e t * Vej)^d 3 vf 0 v 2 cos 2 6 
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where M == - e t M = - e 1 (m/2)/d 3 v f 0 v 2 sin 2 6 is the magnetic moment per unit volume. 
Indicated in parentheses under each term is its ancestry in (16^ >r (18). With a little 
more manipulation we obtain 
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where < v^) is the average drift velocity of guiding centers. As shown by Northrop 

(1961 1, the right hand side of (21) also equals 
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a result which also follows from the moment equation 
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P being the pressure tensor. 
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Although it might have been tempting at first sight of (16) to make off-hand identi- 
cations, such that V B drifts arise from terms containing y f Q , and the curvature drifts 
from those containing (I • yej), this is only partially true in that parts of all these terms 
also lead to the (y X M) x term. Thus when looking at (16). it is difficult to think in 
terms of guiding center motion. 
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III. THE ELECTRIC FIELD IN A RIGIDLY COROTATING MAGNETOSPHERE 


As a preliminary to our using Equation (16) in following papers, we present here a 
short derivation of the form which E must have in order that the magnetospheric plasma 
corotate. The source of B may be finite in size and arbitrary in shape. It is rotating 
about a fixed axis which we take to define the z-axis of a cylindrical “laboratory” 
system, so that H By rigid corotation we mean that all physical quantities are 

time invariant when seen by an observer moving with the velocity fix r, However, if 
the magnetic field is not axisymmetric with respect to U t a 3 B/3t will exist at every fixed 
spatial point £ in the laboratory coordinate system and have a value (Goldstein [19501, 
Backus [1956J, Birmingham and Northrop [1968]) 
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Maxwell’s Curl Equation hence becomes 
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Further 
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where u is the flow velocity of the plasma, assumed to be ideally conducting. 
Next we postulate that 



so that the plasma flow consists of a component of amplitude (pa), to be determined from 
(24), in the “corotational” direction (# may have a field-aligned component) plus the arbitrary 
field-aligned portion k(^). 

Inserting (25) and (26) into (24) and performing the vector differentiations one 
obtains the equation 

d 3B 

— (aB) - a 4 X B - p$ B • Va = SI fl X B. 

dtp - - dip (27) 

Equation (27) is trivially satisfied by a * const. - SI. The “corotational" electric field E 
thus is 


E - - — — | X B (28) 

~ c " 

Note that since k(jr ) in (26) is still arbitrary it could be chosen to be 

k = -pSltp • e j 

so that u is totally perpendicular to B. Our result (27) has the same form as Mestcl’s 
11961 1 solution for the case of a rotating azimuthaliy symmetric magnetosphere whose 
rotation and magnetic symmetry axes are parallel. The form of the expression for E 
thus is unchanged by generalizing to this non-axisymmetric geometry. In applications 
such as in the following paper to Jupiter’s magnetosphere, the electric field has not been 
measured and it is desirable to postulate one in order to reduce the number of unknowns. 
Although (28) undoubtedly is not exact, it is a stage better than assuming the E which 
one would have for a magnetic dipole field. 
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